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ABSTRACT 

We examine the hypothesis that plasma associated with "Type 11" spicules is heated 
to coronal temperatures, and that the upward moving hot plasma constitutes a signif- 
icant mass supply to the solar corona. ID hydro dynamical models including time- 
dependent ionization are brought to bear on the problem. These calculations indicate 
that heating of field-aligned spicule flows should produce significant differential Doppler 
shifts between emission lines formed in the chromosphere, transition region, and corona. 
At present, observational evidence for the computed 60-90 kms~^ differential shifts is 
weak, but the data are limited by difficulties in comparing the proper motion of Type- 
II spicules, with spectral and kinematic properties of associated transition region and 
coronal emission lines. Future observations with the upcoming IRIS instrument should 
clarify if Doppler shifts are consistent with the dynamics modeled here. 

Subject headings: Sun: atmosphere; Sun: chromosphere; Sun: transition region; Sun: 
corona 
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Introduction 



Traditionally, research on the "coronal heating problem" has focused on the transport and 
dissipation of non-radiative energy directly into plasma at heights typical of the solar corona (e.g. 



Kuperus et al. 



1981 



Parkeiill994l : IWalsh and Ireland! l2003l ). 2Mm and more above the visible solar 



surface. At a first glance, this seems odd, since the corona is a tenuous structure fed by mass, 
momentum and energy from below. Indeed, processes operating solely in coronal plasma cannot 
actually produce a corona unless one already exists. Yet five decades of detailed spectral analysis 
from space experiments have implied that there is a significant net downflow of energy in th e form of 
heat conduction and enthalpy from the corona towards the underlying chromosphere (e.g. iMariska 
I992I ). The evidence in support of this picture includes agreement between emission measures, 
above 10^ K, from models dominated by (downward-directed) heat fluxes, and the observation that 
the profiles of transition region lines, while highly variable, are nevertheless mostly red-shifted. In 
addition, much magnetic (free) energy can be stored in the coronal volume in the form of current 
systems and/or waves. Thus the community has become accustomed to working largely within the 
picture of in situ dissipation of mechanical energy directly within the coronal plasma. 



However, observations from the Solar Optical Telescope (SOT, iTsuneta et adl2008l ). mounted 
on the stable, seeing-free platform of the Hinode spacecraft (jKosugi et alhooi ) have revived some 
earlier ideas conc erning the coupling of t he corona to the underlying chromosphere. Data obtained 
and analyzed by Ide Pontieu et al\ (|2007l ) have produced quantitative information about two pop- 
ulations of "spicules", jets of material moving supersonically above the limb and into the corona, 
something that was not ach i evabl e with the previous generation of ground-based observations (e.g. 



Roberts 



1945 



Beckers 



1968 



19721 ) . De Pontieu and colleagues identified a class of finely structured 
spicules with apparent motions far faster ('~ 100 kms~^) than previously measured (~ 20 — 30 
kms~^). These "Type-II" spicules, unlike their longer lived brethren, disappear in a fashion sug- 
gesting that some of the cool material is heated during the spicule's lifetime. Statistical correlations 



ode. TRACE. STEREO, and SDO soacecraft (De Pontieu et al. 


2009: 


Mcintosh and De Pontieu 


2009; 


Mcintosh et al. 


2010: 


De Pontieu et al. 


2011 


)• 



This work suggests that significant energy deposition occurs in the plasma associated with 
spicular events. Recent MHD simulations have provided one scenario in which Lorentz forces and 
Joule heating on a wide rang e of different field lines can explain some of the observed phenomena 



( Martmez-Sykora et al 



2011 



). Here we test the possibility that the heating and acceleration of 
plasma occurs in a field-aligned flow, the freshly heated spicular material filling the overlying 
atmosphere with material heated to higher, perhaps coronal, temperatures. In other words, we test 
whether a simple hydrodynamic approach can reproduce the observed phenomena. These studies 
are of interest because the association of hot plasma with spicules s uggested that "energy d eposition 



at coronal heights cannot be the only source of coronal heating" (jDe Pontieu et al.ll2009l ). 



Ideas along these lines have been around for some time. iThomasI (119481 ) . interpreting Roberts' 
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spicules as hydrodynamic jets, noted that "the directed mechanical energy of the jet be comes con- 
verted to random thermal energy, in part, as the jet moves though the atmosphere". iMivamotc 
(| 19491 ) studied the viscous dissipation of kinetic energy of spicular motions. But with the ob- 
servations the n available, the 30 kms~^ speeds of spicules was insufficient to produce coronal 
temperatures. iPneuman and KoppI (jl978l ) estimated mass fluxes from available data, which with 



mass conservation, prompted them to suggest that "the observed [transition region] downflow repre- 
sents spicular material r eturning to the chromosphere after being heated to coronal temperatures" . 
Athav and Holzeii (|l982l ) noted that "spicular material is raised well above the height that would be 
achieved by a projectile of the same initial velocity, thereby obtaining gravitational potential energy 
much in excess of its initial kinetic energy". Athay and Holzer concluded that "if sufficient heat 
is added to spicules, in conjunction with their acceleration, the spicule phenomenon may also play 
a major role in the production and maintenance of much of the solar corona." In a paper entitled 
"the coronal heating paradox", lAschwanden et al\ ([20071), proposed that the the phrase "coronal 
heating problem" is a misnomer. Instead, they argue that coronal heating should be considered 
as two processes, a "chromospheric heating problem" and "co ronal loop filling proce ss" , moving 
the problem of energy dissipation into chromospheric plasma. iHansteen et al\ (|2010l ) present so- 
phisticated 3D radiation MHD calculations from sub-photosphere to corona and conclude that, for 
calculations driven by the work done on photospheric fields by convection, and dissipated by Ohmic 
heating, "the heating per unit mass [. . . ] necessarily is concentrated toward the transition region 
and low corona" . 



Before proceeding, we clarify some terminology, following lAthay and Holzerl (j 19821 ) . Phrases 
such as "heating in the chromosphere", "heating in the transition region", as well as "heating 
in the corona" abound. We will avoid this usage, as these phrases seem to imply that there is 
something about the thermal properties of these plasmas (e.g. ionization state, conductivities) 
that delineate different mechanisms affecting the three regimes. This may well be the case, but few 
studies deal with such differences (one example that does is the work of iGoodmanI |2004| ) . Also, 
this language suggests, perhaps subliminally, that energy transport mechanisms must be found 
that fit in a pre-existing thermal structure. Instead we see the chromosphere, transition region, 
and corona, including the associated zoo of observed phenomena we call spicules, fibrils, explosive 
events, sprays, blinkers, etc, as observational manifestations of the transport and dissipation of 
mass, momentum and non-radiative energy throughout the Sun's atmosphere. When viewed in a 
more physical sense, one avoids difficulties arising from observationally-based language use, such 
has, how can so much heating appear to be confined to the thin solar transition region? How can 
an observational phenomenon like a spicule can be considered as a source of mechanical heating? 
Why need one consider coronal heating to be paradoxical at all? 
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Calculations 



The work of lDe Pontieu et al\ (j2009l . I2OIII ) suggests that significant energy deposition occurs in 
plasma associated with spicules. Nearly all models on spicules "begin with some form of deposition 
of energy in the photospheric or chr omospheric po rtion of a magnetic flux tube which extends from 
the photosphere into the cor ona" (jSterlingj l200Cl ). a notable exception being the recent work of 
Martfnez-Sykora et al\ (120 111 ). In their detailed calculations, the upward jet motion is essentially 
a pressure-gradient-driven field-aligned flow, initiated by the horizontal compression of plasma by 
the Lorentz force associated with emerging flux causing strong Joule heating in and around the 
jet-like feature. Here, we will study the response of the atmosphere to ad-hoc forcing and heating 
terms. Since the Type II spicules are very thin (high aspect ratio), as a first attempt to compare 
simulated line profiles with observ ations, we perfor m ID gas dynamic calculations solving the set 
of physical equations described bv iHansteenl (|l993l ). These equations follow the time evolution of 
a single fluid's mass, momentum, total energy, and also the number densities of the ground levels 
of atoms and ions of H, He, C, N, O, and Si. The adopted ionization and recombination rate 
coefficients, and the special treatment of ionization of hydrogen, are described in the Appendix. 

Sources and sinks for momentum and energy include gravity and explicit heating per proton and 
neutral hydrogen atom. Radiative losses in the energy equation were treated using the computed 
time-dependent population densities for H, He, C, N, O, Si, losses from other elements and fr ee- free 
emission were treated using a lookup table computed by PGJ using the DIPER package (jJudge 
20071 ). All these c alculations are b ased upon bound-bound collision strengths mostly from the 
CHIANTI project (jLandi et al.ll2006l ). Heat conduction is included using the Spitzer formula. The 
internal energy of the gas is computed using the translational and internal degrees of freedom, i.e. 
the usual perfect gas energy density plus the ionization energy. The internal energy stored only in 
H and He ionization is included. For H, the internal energy is modified by the Balmer continuum 
radiation and so depends on radiative transfer. In the chromosphere we set it to be just | of 
the ionization energy, as j of this energy is supplied by photoionization from the n = 2 levels by 
photospheric radiation. This correction has a minor effect on the calculations. 

The equations, which we solve in their conservative forms, are of mixed character. Advection 
dominates the time evolution of all equations except when the plasma reaches coronal temperatures, 
where the non-linear diffusive heat conduction term domi nates the energy flu x. We therefore 
adopted an operator splitting time integration scheme (e.g. iHansteen et a?.ll2010l ). integrating the 
advection terms using a Lax-Friedrich finite difference scheme, accurate to first order in space and 
time. (We tried slope-limited higher order schemes which are more accurate for smooth solutions. 
But these introduce interpolation errors when, as in the transition region modeled here, gradients 
are very steep). After each advection step, we apply a Crank-Nicolson scheme to integrate up the 
temperature changes due to conduction. A uniform, fixed grid of 4097 points was used (= 10 km 
resolution) along a semi-circular loop of total length 20 Mm (footpoints separated by 40/7r Mm) 
with the dense "photosphere" at both footpoints. Boundary values at the photosphere were held 
fixed. 
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2.1. Initial state 

An initial state was constructed in which protons and hydrogen atoms were heated at the rates 

€p = 2x lO^^^np erg cm~^ s~^ (1) 
ejy = 1 X W^'^nn erg cm~^ s"-^ (2) 

where Up and nn are number densities of protons and neutral hydrogen atoms respectively. These 
values produce a chromosphere and a corona with a peak temperature of lO^K. The equations for 



conser vation of mass, momentum and energy were evolved, starting from model C from lVernazza et al 



(|l98ll ). together with rate equations for the ground states of neutral H and and He, He"*", and for 
H"*" and He^"*". Figure [T] shows the configuration of the starting solution. It is close to equilibrium, 
but has a small yet persistent subsonic downflow in the lower transition region. 



2.2. Forcing 



Guided by the observational analyses of iDe Pontieu et al\ (j2009l . l201ll ) , the initial state was 
forced in an ad-hoc way by accelerating chromospheric plasma and heating it for ti seconds, followed 
by another t2 second period of more intense heating but no acceleration (coasting). Here we discuss 
two calculations, with parameters listed in table [TJ In Calculation "A", plasma initially between 
2000 and 2400 km in height is heated and accelerated for 30 seconds, and then simply heated 
for another 30 seconds. Calculation B has shorter acceleration and heating phases (15 second 
durations), heats/accelerates plasma 100 km lower and higher (roughly one chromospheric scale 
height), but it is heated in the coast phase by a larger factor. 

The heated plasma was initially located at the top of the chromosphere, extending into the 
low corona, and as the calculations evolved the same plasma was heated by tracking the plug of 
plasma. The net upward acceleration of the plug was set to 2 km s~^, about four times the solar 
gravitational acceleration. The heating rates for both neutrals and protons were increased three- 
fold during the first period, to counteract cooling by adiabatic expansion, and then by much larger 
factors during the second period. In this way we mimic the heating of plasma already accelerated to 
~ 35 and ~ 25 kms~^ to coro nal temperatures f or A and B respectively, a little less than speeds of 
typical spicule observations of lDe Pontieu et al\ (j2009l ): iMcIntosh and De Pontieul (|2009l ). In order 
to avoid unphysical discontinuities, the acceleration and heating were smoothed using a Gaussian 
spatial profile of width (FWHM) 330 km about the initial state. 

Figures [2] and [3] show the atmospheric parameters 20 seconds after the beginning of the accel- 
eration, and then 20 seconds after the strong heating began, for calculation A. The sound crossing 
time fo r the calculation shown is abo ut 150s. Note that, unlike the far more sophisticated calcula- 
tion of iMartmez-Sykora et all (j201ll ). these two calculations do not produ ce long, cool spicule- like 
structures, a well known problem with most spicule models (|Sterlingj|2000l ). 
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Our parameter choices parameters are arbitrary in physical terms. But they produce results 
typical of those soluti ons in which chrornosph e ric pl asma is accelerated and heated in the fashion 
broadly suggested by iDe Pontieu et al. (|2oo9l . boiih . Forcing the calculations in denser, deeper 
layers produces warmer chromospheres with no fast outflow to the corona; forcing them higher 
produces insignifcant chromospheric dynamical effects. 



2.3. Results 



First we discuss calculation A. The response of the atmosphere to the input of momentum 
and energy is simple: the plug of forced plasma initially expands mostly outwards into the tenuous 
coronal plasma reaching speeds of ~ 35 km s~^ . After this first phase the additional extreme heating 
provides, through the associated over-pressure, additional outward acceleration of coronal plasma 
to near 90kms~^. In the uppermost regions of the chromosphere, some dow nflow is seen as the hig h 
pressure plug pushes downwards, as discussed in the MHD calculations of lHansteen et al\ (j201Cl ). 
Results for calculation B are similar, except that the extra heating leads to upflow speeds of ~ 120 
kms~^. Also, a cool (~ 3 x lO^K), dense plug of plasma is driven into the corona in calculation B 
behind the shock wave propagating upwards into the corona. 

We have investigated the infiuence of certain assumptions in our calculations. Our results 
are sensitive to the treatment of radiation losses. We treat the losses as effectively thin, as in 
most earlier work, but there are two si gnificant, physically distinc t effects that preclude our use 



of standard radiation loss curves (e.g. iRavmond and SmithI 119771 ). The first is to account for 



photoionization of hydrogen by photospheric Balm er continuum radiation. Non-LTE radiative 
transport is beyond the scope of the present paper (cf. lCarlsson and Steinlll995l ). but the qualitativ e 
effects of Balmer continuum photoionization on radiation losses have been studied by lAthayi (| 19861 ). 
Ionization via Balmer continuum photoionization lowers the radiation losses from hydrogen La, a 
line that otherwise can dominate the plasma radiation losses between 10^ and 2 x lO'^K. The second 
effect is that for the ions treated in detail, we solve radiation losses without assuming ionization 
equilibrium. Our treatment of hydrogen is discussed in the Appendix. The two effects both serve 
to increase the ionization of hydrogen relative to radiation- free optically thin conditions, allowing 
more energy to propagate upwards into the low corona that otherwise would have been radiated 
by the chromosphere. 



2.4. Line profiles 

To compute profiles of selected EUV emission lines, we take the time-dependent populations 
of all ions of H, He, and of all charged ions of C, N, O, and Si^"^ through Si^"^, and use these 
populations to compute the emitted power in selected lines, which are then integrated along a 
given line of sight to compute the line profiles assuming optically thin conditions. The same data 
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are used in the calculation of radiative cooling for these ions. The approach of solving in time 
just for the ionization states of the ions, ignoring excited levels (except for the special case of 
hydrogen), is justified because for the ions of interest, the ion densities are accurately tracked in 
time and space and th e emitted ra diation can be computed post-facto because of the separation 



of timescales involved (jJudgdl2005l ). We compute profiles using Gaussian profiles including both 
thermal and a 10 kms~^ turbulent component to the line widths, and include the Doppler shifts 
of the moving plasma as seen vertically. To make meaningful comparisons with observations we 
average all the emission which, when observed vertically, lies within 2.175Mm (3" as seen from 
earth) of the loop footpoint projected on the solar surface. This length was chose n simply becaus e 



the "transition region" moves up and down significantly in our calculations (cf . iHansteenl Il993l ) , 
and choosing a smaller length scales can produce artifacts as the transition region emission moves in 
and out of this range. (Also, the point spread functions of the best EUV spectrometers would have 
contributions wider than 1"). Figure U] shows the resulting line intensity profiles of a variety of lines 
of He, C, O and Si as a function of Doppler shift and time, for calculation A, in which additional 
heating and acceleration starts at t = 10 sec. In each panel, the element, central wavelength and 
maximum (wavelength-integrated) intensities are listed. The figure shows profiles of lines with 
increasing atomic number, from top to bottom, and with increasing ionization stage- and hence 
typical formation temperature- from left to right. The intensities calculated generally exceed those 
measured in extant low- resolution observations , but given that the filling factor of type II spicules 



may be very small (e.g., Ide Pontieu et al.l 120071 ) . we do not consider this an important difference. 

The overall dynamics is reflected in the lines of Si, from the chromospheric Si II 1526 A line, 
through the Si X 621 A coronal line. The upward acceleration of chromospheric plasma between 
10 and 40s can be clearly traced in the calculations for the Si II line, and in the C II line. It is 
also visible in the coronal lines (C V, C IV, O VII, O VIII, Si X). But in these lines additional 
acceleration due to the pressure gradient in the hotter plasma extends for an additional 30 sec or 
so. 

In fact the chromospheric and coronal lines behave quite differently: the blueshifts of the 
coronal lines are systematically larger than chromospheric and transition region lines by some 50-80 
kms~^ . The line profile of Si II for example shows chromospheric acceleration, but the blue-shifts 
disappear with the onset of the additional heating at t = 40,s. In contrast, the coronal lines are 
blueshifted relative to the Si II line when t > 40s. After the forcing period the acoustic wave 
generated propagates away from the footpoint, and all the lines become dimmer. In addition, all 
chromospheric lines return near t = 100s to zero Doppler shift and even become red-shifted (faintly 
seen in the lines of helium, for example). 

The different Doppler shifts calculated for chromospheric and coronal lines, found in this 
particular calculation, is a more general result. It is seen in calculation B, where the stronger 
heating rates produce sudden pressure-driven acceleration, so that the coronal line blueshifts exceed 
the Doppler shifts of chromospheric lines by ~ 100 kms~^ (Figure [5]). Several other calculations 
with different heating rates higher and lower in the chromosphere produce qualitatively similar 
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results. 



2.5. Other spectral signatures of the dynamics of heated outflowing plasma 



Another systematic difference is seen in the emergent spectra we have calculated. The spec- 
tral lines we selected include some atomic transitions which have sensitivity to ionization non- 
equilibrium effects, through their relatively high excitation energy. Excitation of helium atoms 
and ions requires a large energy (compared with typical thermal energies found under ionization 
equilibrium conditions) because the excited levels have a larger principal quantum number n than 
the ground levels. Any excitation requires a "An > 1" transition. In contrast, other atomic ions 
can have both An = and An > 1 transitions. The plotted 312 and 457 A lines of C IV and Si IV 
respectively are An = 1 transitions (2s — 3p and 3s — 4p), unlike the 1548 and 1393 A lines which 
are 2s — 2p and 3s — 3p transitions respectively. 

In Figure[6]we show the profiles computed with the full dynamics, but assuming instantaneous 
ionization equilibrium (IE). Comparison with Figure [H clearly reveals departures from IE. The 
IE calculations are usually more variable as it takes some time for ions to become ionized or to 
recombine in response to changes in temperature and density. But, systematically, the An = 
1 transitions are all brighter than IE would predict. This effect results from plasma becoming 
ionized, under conditions where excitation and radiative decay of the levels occurs before the ion to 
which they belong becomes sig nificantly ion i zed. In helium lines this has been invoked to explain 

Pietarila and Judgell2004l : 



their anomalous brightness (e.g. 



Jordan 



1975 



Laming and Feldman 



1992 



Judge and Pietarilal l2004l : I Judgd l2005l ) . A similar effect is seen in the C IV and Si IV ions, but 
here we can compare directly the intensities An = 1 with the An = transitions. In both cases, 
when the emitting plasma is undergoing heating, the An = 1 transitions become enhanced relative 
to the An = transitions by a factor of 3.1 and 2.5 for Li- like C IV and Na-like Si IV respectively. 

Such large departures from ionization equilibrium should be measurable. Simultaneous obser- 
vations of the An = 0, 1 transitions in C IV and Si IV have presumably been obtained using the 
CDS and SUMER instruments on the SOHO spacecraft, but it is beyond the scope of the present 
paper to identify and analyze such data. 



3. Discussion 

Our ID hydrodynamic calculations attempt to capture the evolution of the emitted spec- 
trum of a plug of accelerated and heated plasma flowing along magnetic field lines extending 
upwards from the upper chromosphere. The dynamics and emergent spectra are computed using 
full non-equilibrium in the atomic rate equations. The calculations reveal that the coronal emission 
which accompanies the heating needed to produce coronal plasma from the cool spicular plasma is 
blueshifted by 50 — 90 kms~^ relative to the chromospheric emission (compare lines of He I, C II 
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Si II with the coronal hnes in figs. This is because, in the hydrodynamic regime, the heating 
produces an over-pressure and a subsequent expansion of the plasma into the corona. Thus, if 
indeed the chromospheric plasma in Type II spicules is injected into the corona, at the same time 
being heated in this fashion, then we must expect the velocity distributions of chromospheric and 
coronal lines to differ. Accompanying the prediction of Doppler shift differences, changes in the 
intensities of line ratios that are traditionally viewed as being sensitive to electron temperature, 
are also predicted in our rapidly heated outflow models. These intensity changes result from the 
slowness of ionization/recombination processes relative to dynamical timescales. It would seem 
worthwhile to investigate at omic systems of the Li and Na-like isoelectronic sequences observation- 
ally, following earlier wo rk (IHeroux aLlll972l ') as well as other atomic systems with temperature 



sensitive line ratios (e.g. iPinfield et at 



199 



Such a systematic difference between cool and hot flows associated with Type II spicules has not 
been reported in the observational literature, but we note that the observations have significant 
uncertainties. The c ool components of Typ e II spicules at the limb show apparent velocities of 
order 50-100 kms~^ (Ide Pontieu et aLl 120071 ). with rapid blues hifted events, thought be their disk 



counterparts, showing line-of-sight velocities of order 50 kms ^ (Rouppe van der Voort et aLll2009l ) 



and apparent velocities of order 75 kms~^ (jPe Pontieu et al.ll201ll ). The velocities of coronal 
counterparts are also difficult to determine precisely. They are derived from weak blueshifted 
emission components, seen at the base of much stronger emission cores, near the limit at which EUV 
spectrographs can operate (e.g., instrumental broadening, signal-to- noise, ...). These measurements 
suggest line-of-sight velocities of order 50-150 kms^'^ for the coronal counterparts of spicules. 
Because both the chromospheric and coronal measurements are significantly impacted by viewing 
geometry, radiative transfer and instrumental limitations, it is not yet clear that the observed 
velocity distributions are sufficiently well determined to reveal a velocity difference of 50 kms~^ 
between the chromosphere and coronal counterparts of spicules. 

If the large relative shift in the Doppler motions predicted here are not confirmed by future 
high resolution measurem ents (such as those from the IRIS missiorjl l) , we would conclude that the 
suggestion put forward by Ide Pontieu et all (j2007l ): lDe Pontieu et all (j2009l . l201ll ). namely the close 
correspondence of velocities from spicules with those of Doppler shifted lines in the transition region 
and corona, cannot credibly be accounted for in the type of ID hydrodynamic flow modeled here. 
One possible resolution might simply be that the magnetic fields within the tube are significantly 
twisted such that the spicular plasma has an azimuthal as well as axial component as the Lorentz 
force directs the plasma along field lines. In this case the Doppler shifts computed should be 
multiplied by the cosine of the magnetic pitch angle. 



The recent MHD calculations of lMartmez-Sykora et al\ (|201ll ) might help resolve the problems 
identified here, in that the system modeled is far more complex than our simple calculations and 
does not seem to show such large velocity differences between chromospheric and coronal lines. In 



^http: / /iris. lmsal.com/index. htm 
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their work, they find that chromospheric material is injected into the corona via, first, compression 
due to horizontal Lorentz forces associated with emerging flux, and second, via the field-aligned 
gas pressure gradient. As the field-aligned flow progresses it is heated via Joule dissipation, in a 
fashion not dissimilar from the ad-hoc calculations we present. However, in the modeled MHD 
system, the acceleration and heating acts, with varying strength, on plasma that occurs on an 
ensemble of neighboring field lines. As a result, the coronal emission has contributions from many 
magnetic field lines with varying mixes of chromospheric and coronal plasma/flows. Thus, the line 
profiles of coronal plasma must be carefully synthesized to find the relative contributions of these 
different components. Notably though, the cool plasma that actually flows into the corona in these 
MHD calculations is indeed a field-aligned flow. Thus, if our parametrized heating and acceleration 
profiles are representative of what occurs in the complex 3D configuration, the heated plasma in 
that part of the calculated flow would be expected to behave like the calculations presented in the 
present paper. Current EUV observations lack the spatio-tem poral resolution to reso l ve the finely 
structured dynamics and energetics predicted by the model of iMartmez-Svkora et all (j201ll ). 



The potential problern of Doppler shift differences may be avoided in an entirely different 



scenario (jJudge et al. 



2011 



where it was suggested that Type II spicules, in particular those with 
apparent upward velocities in excess of 50 kms~^, correspond not to jets but to the line-of-sight 
superposition of warped current sheets modulated by AlfvenicH fluctuations driven from below. 
Low-frequency Alfyenic -type wave motions have been reported before in movies of Type II spicules 
(|De Pontieu et aLl 120071 ). Should Type II spicules be proven to be sheets, the calculations here are 
largely irrelevant, and the physical connection between the chromospheric, transitior i region and 



coron al emission proposed earlier would be called into question. But, as suggested by I Judge et al 



(j201ll ) the ID flow picture might apply to a subset of observed spicules with plasma genuinely 
flowing at lower speeds. However, several questions regarding the overall interpretation of spicule 
properties at the limb, their relation to RBEs, and statistical analyses would then arise. 

Whatever the out come, we conclude tha t the role of Type II spicules in supplying mass and 



energy to the corona (jde Pontieu et al\ 120071 : iDe Pontieu et al.ll2009l . l201ll ). is a subject ripe for 



future study. Further observational work on spicule properties and their relationship to the prop- 
erties of the associated hotter plasma seems warranted. The upcoming IRIS mission seems well 
positioned to address this question. 



We are grateful to Paul Cally and an anonymous referee for a careful reading of the paper and 
for helpful comments. 



"^Simply meaning motions dominated by the balance between magnetic tension and inertial terms. 
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A. Atomic rate equations 



The system of equations integrated in time includes equations for the number densities of 
atoms and atomic ions. In ID, these take the form 



N 



n. 



N 



(Al) 



where t, x are time and distance respectively, rii is the population density of atomic state i, and the 
coefficients Pji represent the transition r ate, units s~ ^ , from state j to state i (and vice versa) . Only 
long-lived states need tracking in time (j Judgd l2005l ) . so that the only states we solve for are the 
ground states of atomic ions and bare atomic nuclei. Given the solar radiation field and local thermal 
properties in the corona, the largest contributions to coefficients P are ionization and recombination 
coefficients, by electron impact (charge transfer collisions with H and He a,re sometimes important 
but ar e not included here) . Here we adopt the ionization rate coef ficients oflArnaud and Rothenfiug 



(1985), and fits to recombination rates in the form tabulated by ShuU and Steenbergl ( IQsd ). com- 
puted by one of us (PGJ) from the detailed photoionization cross sections computed by the OPAC- 
ITY project (ISeatonlll987l). These differ su bstantially from, and should be more accurate than, the 
calculations of (jshull and Steenberg 1982 ). A com parison of recombina t ion r ates compute d from 
the O PACITY project and more detailed work of iNahar and PradhanI (j 19921 ) is given in (jjudge 
20071 . Fig. 13jf|. The systematic differences show the rates to agree within 0.2dex with some 
excursions of a factor of 2 possible. 

Hydrogen ionization occurs not quite so simply as described above, and it is treated differently. 
This is because the n = 2 levels can have a significant population owing to scattering in Lya within 
the chrom osphere, and Balmer continuum radiation from the up per photosphere can photoionize 



hydrogen (|Hartmann and MacGregorlll980l : IVernazza et a?.lll981l . e.g.). Thus we include photoion- 
ization from the n = 2 levels. We compute the n = 2 populations assuming that Lya is in detailed 
balance in the chromosphere (mostly neutral) but that the population is greatly reduced below this 
limit in the corona. As an ad-hoc, qualitative representation of these effects we apply the non-LTE 
factor 



n2 
ni 



ni 



(A2) 



0-^ ym + n^^ 

where the asterisk refers to LTE populations, and is the p opulation density of protons. We 

We assume 



adopt a photoionization rate of 8 x 10 s from the n = 2 level (jVernazza et al. 



1981 



radiative detailed balance throughout in the Lyman continuum and thus set photoionization and 
recombination rates to zero in that transition. This is not a serious problem in that it is correct in 
the chromosphere, and photoionization is not very important throughout the transition region and 
corona, relative to electron impact ionization. 



^Available at |http://www.hao.ucar.edu/modeling/haos-diper/[ 
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Fig. 1. — The initial, near-equilibrium state of the atmosphere is shown. 



This preprint was prepared with the AAS IM^^X macros v5.2. 
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Fig. 2. — The state of the atmosphere after 20 seconds of acceleration and moderate extra heating 
in the upper chromosphere is shown, for calculation A. 



Table 1. Parameters of the acceleration and heating 



Calculation 


h 


t2 




Z2 


hi 


h2 




s 


s 


km 


km 






A 


30 


30 


2000 


2400 


3 


30 


B 


15 


15 


1900 


2500 


3 


100 



Note. — ti is the duration of the acceleration and moderate heating phase, t2 the coasting and 
large heating phase. The initial atmospheric plasma between heights zi and Z2 is the region subject 
to the acceleration and heating, hi and /12 measure the factors by which the nominal heating rates 
(equation [1]) are increased during ti and t2 respectively. 
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Fig. 3. — The state of the atmosphere 20 seconds into the intense heating phase is shown, for 
calculation A. 
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Fig. 4. — Spatially-averaged emission profiles from near the accelerated and heated footpoint of 
the loop, for calculation A. The emission was integrated along the vertical and averaged over all 
spatial positions within 2.175 Mm (corresponding to a 3 " spatial bin) of the footpoint in the low 
chromosphere. The panels are labeled with the atomic ion, central wavelength, and computed 
maximum total intensity ( ergs cm~^ s^^ sr^^). The horizontal dotted lines demark the periods of 
acceleration and heating (10 to 40s), and of more intense heating (40 to 70s). 
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Fig. 5. — Same as Figured! but for calculation B, which has higher heating rates in denser plasma 
than calculation A, and shorter periods of heating and acceleration. 
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Fig. 6. — Same as Figure [H but for calculations assuming instantaneous ionization equilibrium. 



